Particle Physics on Ice: 
Constraints on Neutrino Interactions Far Above the Weak Scale 
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Ultra-high energy cosmic rays and neutrinos probe energies far above the weak scale. Their 
usefulness might appear to be limited by astrophysical uncertainties; however, by simultaneously 
considering up- and down-going events, one may disentangle particle physics from astrophysics. 
We show that present data from the AMANDA experiment in the South Pole ice already imply an 
upper bound on neutrino cross sections at energy scales that will likely never be probed at man-made 
accelerators. The existing data also place an upper limit on the neutrino flux valid for any neutrino 
cross section. In the future, similar analyses of IceCube data will constrain neutrino properties and 
fluxes at the 0(10%) level. 
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Ultra-high energy cosmic rays have been observed with 
energies above 10 10 GeV, implying collisions with terres- 
trial protons at center-of-mass energies above 100 TeV. 
Ultra-high energy cosmic neutrinos, so far undetected, 
are expected to accompany these cosmic rays. These 
cosmic neutrinos are especially interesting, because their 
known interactions are so weak that they are highly sen- 
sitive to new interactions at energies far above the weak 
scale, where new interactions are expected in many ex- 
tensions of the standard model (SM) of particle physics. 
In addition, ultra-high energy neutrinos are unique mes- 
sengers, as they are expected to escape from (and point 
back to) even the most dense astrophysical sources. 

The promise of ultra-high energy neutrinos might ap- 
pear to be severely limited by astrophysical uncertainties. 
Event rates constrain only a combination of fluxes and 
cross sections, and so astrophysical uncertainties cloud 
particle physics implications and vice versa. However, 
the event rates for up- and down-going neutrinos depend 
differently on neutrino cross sections |l|,|2(. By combining 
both up- and down-going data one may therefore disen- 
tangle particle physics from astrophysics and constrain 
both the properties of astrophysical sources and the in- 
teractions of neutrinos far above the weak scale. 

Here we consider neutrino telescopes operating under 
ice at the South Pole @, Q . We show that current data 
from the AMANDA South Pole telescope [j| already sig- 
nificantly constrain neutrino cross sections at center-of- 
mass energies yfs ~ 6 TeV, and future data will provide 
0(10%) determinations. These results will be complc- 
mented J|j by future data from the Pierre Auger Obser- 
vatory pfl . The energies probed are far above the HERA 
domain -Js ~ 500 GeV, the highest accelerator energy at 
which even indirect tests of neutrino cross sections have 
been made. 

Simple parton model predictions for neutrino-nucleon 
cross sections :Sj may be suppressed, for example, by 
saturation effects at small x |9|. Such effects have been 
proposed to slow down the power law scaling of cross 
section with neutrino energy to comply with the Frois- 
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sart bound [lfj. On the other hand, neutrino cross sec- 
tions may also be enhanced, for example, by the exchange 
of towers of Kaluza-Klein gravitons [llj , black hole pro- 
duction 0,0, TeV-scale string excitations [l3| . or elec- 
troweak instanton processes 14]. Our results constrain 
all of these possibilities. 

The results derived below also constrain the extrater- 
restrial neutrino flux, which is at present unknown. Neu- 
trino sources may be conveniently characterized as either 
optically thin or thick. In optically thin sources, the nu- 
cleons responsible for neutrino production (through pho- 
toproduction or pp collisions) escape the source and con- 
stitute the observed cosmic rays. The observed cosmic 
ray and diffuse neutrino fluxes are then correlated |l5j |. 
In contrast, in optically thick sources, the neutrino pro- 
genitors do not escape, thus vitiatin g th e relation be- 
tween cosmic ray and neutrino fluxes |l6j . In principle, 
this permits a very large enhancement to the neutrino 
flux. We derive upper bounds on this enhancement that 
are valid for any neutrino cross section. 

We will derive bounds without assuming particular 
neutrino fluxes or cross sections. It will be convenient, 
however, to present results relative to standard reference 
values. For the reference cross section, we choose ctsm, 
the charged current (CC) neutrino-nucleon cross section 
of the simple parton model [8j . For the reference flux, we 
adopt the Waxman-Bahcall (WB) flux <^ B 15]. This 
flux is that of an optically thin source. At production, 
the WB flux has flavor ratios : v e : v T = 2 : 1 : 0, but 
this quickly transforms to 1 : 1 : 1 through neutrino os- 
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for each flavor. In what follows, we focus on neutrino en- 
ergies in the range 10 7 GeV to 10 75 GeV, where the back- 
ground from atmospheric neutrinos is negligible, but the 
extraterrestrial flux is expected to be significant. Here- 
after, we take (E) = 10 7 ' 25 GeV. 

Neutrinos are detected in neutrino telescopes when 
they create charged leptons or showers through CC or 
neutral current (NC) events near the neutrino telescope, 
and the resulting leptons or showers propagate into 
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the experiment's instrumented volume. For down-going 
events, the probability of neutrino conversion is always 
small, barring extraordinary enhancements to neutrino 
cross sections. Letting <7„at denote the total (NC + CC) 
neutrino-nucleon cross section, the down-going event rate 
is therefore 



own 



= c 



down 



0~vN 



(1) 



where <fi v indicates the average extraterrestrial v + T> flux 
per flavor in the energy bin of interest. The constant 
Cdown depends on exposure and acceptance and varies 
according to neutrino flavor from experiment to experi- 
ment. 

For up-going events, the dependence on cross section 
is completely different. At energies above 10 6 GeV, neu- 
trino interaction lengths become smaller than the radius 
of the Earth. For v e and v^, this implies that most 
upward-going neutrinos are shadowed by the Earth, and 
only those that are Earth-skimming |17| , traveling at low 
angles along chords with lengths of order their interac- 
tion length, can produce a visible signal. For v T , this 
attenuation is somewhat offset by effects of regeneration 
► t ► z/j- jl8[ . However, in the energy range of in- 
terest, the expected events will be largely dominated by 
neutrinos surviving until the region of the detector. 

The dependence of upward-going event rates on 
anomalous neutrino cross sections depends on the source 
of the anomaly. We consider two prominent cases. First, 
in many new physics cases, the SM CC cross section <tsm 
is not altered, but there are new neutrino interactions 
that produce showers. Letting o v n be the total neutrino 
cross section, including the standard model CC and new 
physics contributions, the up-going event rate is |(J 
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As discussed in Ref. gj, Eq. J5J) is valid assuming L l <C 
V < i?^, where L l is the typical lepton path length in 
Earth, V is the neutrino interaction length in Earth, 
and R® is the radius of the Earth. This corresponds 
to E > fO 7 GeV. Here, <j u n indicates any enhancement 
of the cross section which will increase the event rate for 
down-going neutrinos, but because of absorption will sup- 
press the up-coming events. The latter can be achieved 
through cuts on shower energy fraction greater than or 
equal to that characterizing the CC SM process. Extreme 
enhancements to <7„jv may reduce L v to L , leading to a 
different parametric dependence in Eq. J2J, but we ne- 
glect such cases here. Second, we consider the possibility 
of screening, in which both the standard model NC and 
CC cross sections are reduced equally. In this case 0,0, 
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where <j v n an d ^sm are CC cross sections with and with- 
out screening, respectively. 



Given the parametric dependences of Eqs. QJ, 10), and 
(|3*|). we now consider the implications of existing data 
from AMANDA. Assuming the standard model CC neu- 
trino interaction, the AMANDA Collaboration has de- 
rived the 90% CL upper bound on the diffuse neutrino 



flux E 2 < lax = 3.3 x 10 



7 GeV cm 2 s 1 sr 1 per fla- 



vor 0], assuming an E~ 2 dependence of the flux, valid 
for 10 6 GeV to I0 9 5 GeV neutrinos. Since the energy 
distribution of the AMANDA data peaks in the energy 
bin of interest, it is reasonable to use as the 

upper limit in the bin. Here we generalize this bound to 
the case in which there are two unknown quantities: the 
neutrino cross section and the neutrino flux. Applied to 
v e down-going events, the constraint implies 
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Dividing Eq. Q by gi ves 



< 16 



(4) 



(5) 



V\¥B CT SM 

at 90% CL. A similar analysis for up-going events yields 
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for the case of new physics contributions, and 
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for the case of screening. 

These constraints exclude the shaded regions of Fig. ^ 
The upper region is excluded by down-going data and 
the lower region is excluded by the up-going data, as- 
suming screening. The upper and lower regions meet at 
<3i>n /^sm = 1. As a result, for any neutrino cross section, 
we find an upper bound on the neutrino flux in the en- 
ergy range 10 7 GeV to I0 7 5 GeV of <j)" < 16^ B . Note 
that the lower shaded region limits the amount by which 
screening effects can suppress ctsm- 

How will these results improve in the near future? We 
now consider the possible implications of IceCube, the 
successor experiment to AMANDA. In the energy range 
of interest the r decay length is comparable to the instru- 
mented volume; thus, one can observe all v r topologies: 
a t track followed by the r-decay shower ("lollipop"), 
a hadronic shower followed by a t track which leaves 
the detector ( "popillol" ) , and double bang events. All 
of these distinctive topologies allow a direct and precise 
measurement of the incoming neutrino energy |19j . We 
also note that the absence of oscillation precludes a v T 
atmospheric background. For these reasons, the detec- 
tion prospects for up-going neutrinos are brighter for v T 
than the other flavors, and we focus on them below. 

To evaluate the prospects for IceCube, we must deter- 
mine Cdown and C up for IceCube. We focus on the case 
in which new neutrino interactions modify the NC cross 
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sections, but leave the CC cross sections invariant. In 
NC processes most of the energy is carried off by pions. 
At TeV energies, the interaction mean free path of tt^ 1 
in ice is orders of magnitude smaller than the pion decay 
length, and so nearly all energy is channeled into elec- 
tromagnetic modes through ir° decay. To estimate the 
efficiencies for down-going events in the NC channel, we 
therefore adopt as our basis of comparison the electro- 
magnetic showers induced by v e | 20] . With this in mind, 

Cdow„ = 2ir(j > ^ B <j SM N T TAE f 7?(cos0) dcosO , (8) 

Jo 

where iVx is the number of target nucleons in the effective 
volume, 0.5 < rj(cos9) < 1 is the experimental efficiency 
for detection of electron neutrino showers with zenith an- 
gle 9, T is the running time, and AE = 2.2 x 1 7 GeV. 

IceCube is now under construction, in the process of 
growing to its final size. Including absorption effects, 
it will have an effective aperture (ACl) e g ~ tt/2 km 2 sr 
for detecting v T pj. The detector will consist of 80 
kilometer-length strings, each instrumented with 60 op- 
tical modules. The number of target nucleons may 
therefore be estimated to be Nt ~ 6 x 10 38 . To re- 
main conservative, we take J Q T](cos6) dcosO = 0.8. In- 
serting these numbers into Eq. (JBJ), along with a life- 
time of the experiment T — 15 yr and ctsm ((E)) — 
2 x I0~ 33 cm 2 0! we obtain C doW n ^ 4. The cor- 
responding quantity for up-going events in our energy 
interval is C up = (AO) eff T f (ff (E)dE ~ 20, where 
J (j) T {E)dE = 8.5 x lO" 1 km~ 2 yr" 1 sr" 1 is the r-lepton 
flux emerging from the Earth due to (unabsorbed) v r in- 
teractions inside the Earth for incoming <jf = 4>wb an d 
&vN = csm EH- 

Given these estimates of Cdown and C up , we now deter- 
mine projected sensitivities of IceCube to neutrino fluxes 
and cross sections. For a given set of observed rates 
7V°p S and A/Jo Wn , two curves are obtained in the two- 
dimensional parameter space by setting A/"°p S = A/ up and 
•^down = Miown- These curves intersect at a point, yield- 
ing the most probable values of flux and cross section 
for the given observations. Fluctuations about this point 
define contours of constant % 2 in an approximation to a 
multi-Poisson likelihood analysis |22| . 

In Fig. ^ we show results for two representative cases 
that are consistent with the AMANDA bounds derived 
above. In the first case, we assume <t v n — ^sm and 
cjf = 0wB' leading to 4 down-going and 20 up-going 
events. The 90%, 99%, and 99.9% CL contours are those 
given in the lower left of the figure. (These contours will 
be slightly distorted for <7 u n/&sm < 1, where Eq. (J2J 
receives corrections, but we neglect this effect.) We see 
that, even in the case that event rates are in accord with 
standard assumptions, the neutrino-nucleon cross section 
is bounded to be within 40% of the SM prediction at 
90% CL. This is at a center-of-mass energy y/s ~ 6 TeV, 
far beyond the reach of any future man-made accelera- 
tor [23. In the second case, we consider a scenario in 




FIG. 1: Projected determinations of neutrino fluxes and cross 
sections at y/s ~ 6 TeV from future IceCube data. Two 
cases with (A/J^ n , A/j£ s ) = (4,20) and (35,20) are consid- 
ered. In both cases, the best fit flux and cross sections are 
shown, along with the 90%, 99%, and 99.9% CL inclusion con- 
tours. Contours of constant 7V up = 20 (dotted), A/down = 4 
(left dashed), and A/down = 35 (right dashed) are also shown. 
Neutrino fluxes and cross sections excluded by AMANDA at 
90%CL are also indicated: The upper (lower) shaded region is 
excluded by null results for down-going (up-going) events [f|- 



which the number of observed upcoming events remains 
at 20, but the number of down-going events is 35. In the 
second case, clearly one has discovered new physics at 
well beyond 5a. 

It is noteworthy that in the energy bin of interest, 
the predicted [24[ diffuse flux of neutrinos from a uni- 
form distribution of blazars, <^ GN (iJ), is about 9 times 
larger than the WB flux, i.e., ^agn((^)) — 9^wb- These 
neutrinos are expected to be produced in optically thick 
cores of blazars when ultra-high energy protons scatter 
off the accretion disk orbiting the active galactic nucleus 
(AGN) [25|. By re-scaling the integrated luminosity in 
Eqs. ||TJ and 0, it is straightforward to see that if the 
extraterrestrial flux is at the 4>agn ^ eve K i> n V r °f oper- 
ation IceCube will probe 40% (70%) enhancements from 
SM predictions at the 90% (99.9%) CL. 

Finally, we note that the AMANDA constraints have 
significant consequences for what may be seen at Ice- 
Cube. Substituting Eq. ||SJ| into Eq. (Q) with Cdown de- 
termined for IceCube, the v e down-going event rate in 
the energy range 1 7 GeV to 10 7 - 5 GeV is constrained to 
A/"down < 4 yr" 1 at 90% CL .26j. Event rates at IceCube 
for low scale gravity models have been presented p^ . 
These event rates, along with this bound on A/down, can 
be used to constrain the multidimensional Planck scale. 

In summary, we have shown that the sensitivity of 
neutrino telescopes in the Antarctic ice to both up- 
and down-going ultra-high energy neutrinos provides 
a powerful probe of ultra-high energy neutrino fluxes 
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and anomalous neutrino interactions. Current null re- 
sults from AMANDA already provide interesting con- 
straints on the flux-cross section parameter space. First, 
these results constrain both suppressions of the neutrino 
flux from putative screening effects and enhancements 
from new physics. Second, they exclude large fluxes 
(j) v > 160^ B for neutrino energies between 10 7 GeV and 
10 7 ' 5 GeV, for any neutrino cross section. These energies 
correspond to neutrino-nucleon collisions at y/s ~ 6 TeV, 
far above the weak scale and likely never to be accessible 
at man-made accelerators. 

In the future, IceCube will be able to determine both 
neutrino fluxes and cross sections with impressive accu- 
racy. We have considered neutrinos with energies be- 
tween 10 7 GeV and 10 7 5 GeV. For standard model 
cross sections and the WB flux, 40% (70%) enhance- 
ments from standard model predictions may be excluded 



at 90% (99.9%) CL in 15 years of running. Should the 
neutrino flux be at the level predicted for optically thick 
blazars, these bounds can be attained after two years 
of data collection at IceCube. Our analysis assumes an 
extraterrestrial neutrino flux with 1:1:1 flavor com- 
position. If IceCube finds a different flavor mix (2^], it 
will not be difficult to repeat this analysis for the correct 
flavor ratios. 
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